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We study the structure of the time evolution of the density matrix in contact with a 
thermal bath in a standard projection operator scheme. In general, the equation of motions 
of the density matrix with dissipative effects tends to lead any initial state into a steady 
state. The reduced density matrix of the system in the steady state is obtained by tracing 
out the degree of freedom of the thermal bath from the density matrix for the equilibrium 
state of the total system. This reduced matrix is modified by the interaction, and is different 
from that of the equilibrium of the system alone. In a commonly used equation of motion 
of the reduced density matrix, we have three terms, i.e., a term of the quantum mechanical 
evolution of the system density matrix, a term of the non-Markov evolution due to memory 
effects, and a term depending on the initial total density matrix. We make clear roles of 
the three terms by explicit calculations of the contributions of the terms to the steady state 
density matrix. By making use of the role of each term, the properties of the commonly used 
quantum master equation are examined. For example, if we do not include the imaginary part 
of the second term, the quantum master equation leads the reduced density matrix into that 
of the equilibrium of the system without modification. On the other hand, if we include the 
imaginary part, the steady state density matrix of the system satisfies the master equation 
up to the second order of the interaction. This property indicates that, in the representation 
that diagonalizes the system Hamiltonian, the leading order (the second order) terms of the 
off-diagonal elements of the steady state solution of the master equation agree with those of 
the density matrix of the equilibrium state. 

KEYWORDS: quantum master equation, projection operator, coarse-grained dynamics, relax- 
ation phenomena 

1. Introduction 

Recently, it has become important to understand explicit quantum dynamics in many 
fields. For example, in the so-called single molecular magnets, the quantum dynamics of mag- 
netization has been extensively studied. 1-10 In these systems, the energy levels are discrete, 
and the quantum mechanical dynamics of a wave function under a sweep of external field is 
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described by the combination of adiabatic and non-adiabatic transitions, where the Landau- 
Zener formula plays an important role. In addition to the pure quantum mechanical dynamic, 
effects due to contact with a thermal bath also give important contributions. It has been 
pointed out that the magnetization shows a plateau in the magnetization process under a 
swept magnetic field. 6 ' 7 ' 11 To understand this phenomenon, a kind of rate equation for the 
populations of energy levels 6 and a quantum master equation 7 have been studied. It turned 
out that the appearance of the plateau is a generic phenomenon in systems in which the 
state changes almost adiabatically with a small inflow of heat from a thermal bath. Quantum 
master equation approaches have also been used for various properties of the magnetization 
processes of single molecular magnets. 8,9 Recently, the coherent driven Rabi oscillation has 
been observed in V15, 10 where decoherence processes play important roles. The importance 
of quantum mechanical dynamics has also been pointed out in electron motion under a time- 
dependent field, 12 and in various operations related to quantum computing and others. 13-16 
In this system, dissipation due to contact with the environments also play important roles. 

To describe relaxation phenomena in quantum dynamics, we may derive the equation of 
motion of the density matrix of the system using a method of the projection operator. 17-37 
In general, we derive the equation from the Bloch equation of the total system. The total 
system consists of the system (S), thermal bath (B), and interaction between them (I), the 
Hamiltonians of which are given by Tis, Hi, and Tt-Q, respectively: 

TCt = TLs + XH\ + Ti.B, (1) 

where A is a constant that controls the strength of the interaction. The dynamics of the density 
matrix of the total system is given by the Bloch equation 

ih^ = [H T ,pr}. (2) 

Generally, we assume that the contact with the thermal bath causes dissipation effects, and 
that the system would relax to the thermal equilibrium state. 36,3 To describe such processes, 
we adopt an equation for the time evolution of a reduced density matrix in which the degrees 
of freedom of the thermal bath are traced out. In general, the equation leads an initial state 
into a steady state. This steady state should be compatible with the equilibrium of the total 
system. The reduced density matrix of the system in the steady state is obtained from that of 
the total system by decimating the degrees of freedom of the thermal bath, and it must be the 
steady solution of the time evolution of the reduced density matrix. Because of the interaction 
between the system and the thermal bath, the reduced density matrix of the steady state is 
modified from that of the equilibrium of the system itself. 

The equation obtained by the projection method consists of three terms: a term of the 
quantum dynamics due to the system Hamiltonian, a term of the non-Markov evolution due 
to memory effects, and a term depending on the initial total density matrix. We study the 
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contributions of these three terms to the steady state explicitly. We point out that the third 
term plays an intrinsic role in the steady state, and that artificial choices for the third term, 
e.g., a product form of the system part and thermal bath part, would not make the steady 
state simpler. 

Moreover, we discuss the properties of the quantum master equation in the context of the 
equation of motion of the reduced density matrix. A conventional quantum master equation 
leads the reduced density matrix into that of the equilibrium of the system itself, while the 
time evolution of the reduced density matrix leads the reduced density matrix into a modified 
one owing to the interaction with the thermal bath. The modified one agrees with the density 
matrix of the system interacting with the bath in the leading order of interaction. Namely, 
the second order terms of the off-diagonal elements of the solution in the representation that 
diagonalizes the system Hamiltonian correctly reproduces those of the equilibrium density 
matrix of the system interacting with the bath. We study the roles of the real and imaginary 
terms of the second term of the equation of motion of the reduced density matrix, and clarified 
the difference in the steady state. 

2. Formulation 

2. 1 Equation of the density matrix 

First, let us briefly review a standard derivation of an equation for the reduced density 
matrix of the system p$ by tracing out the degree of freedom of the bath from the density 
matrix of the total system px ;17 

PS = Tr B/ o T , (3) 

where Ttb denotes the trace operation for the degree of freedom of the thermal bath. Here, 
we adopt a projection operator V to the Hilbert space of the system as 

Vp T = PbTtbPt, (4) 

and 

V' = 1 - V, (5) 

where pb is the equilibrium density matrix of the thermal bath. The equation of motion of 
the reduced density matrix VpT is given by the following three parts: 

= ViCVpT + ViC f e^ v ' lC V'iCVp T {r)dT + ViCe^-W^V p T (t ), (6) 

where 

i£ PT = — [H s + XHi + H b ,Pt}. (7) 

This equation consists of three terms: the first term denotes the quantum dynamics due to the 
system Hamiltonian, the second term represents the non-Markov evolution due to memory 
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effects, and the last term depends on the initial total density matrix pr(to). Here it should be 
noted that this equation is part of the exact equations for the total density matrix together 
with the equation of V'pt- Thus, we cannot regard this relation as an equation of motion of 
ps, but a relation concerning p$ in terms of /5t(*o)- 

Here, we study cases with a weak coupling, and consider the dynamics up to the second 
order of A. Then eq. © is expressed in the form (see Appendix A) 

dps _ J_ 
dt ih 



n s ,ps] +Tr B iA f e {t - T)iCa i£ipBPs(r)dr 
Jt 



+ATr B i/:ie (< " <o)i£o P / pT(io)+A 2 TrBii:ie ( ^* o)4£o f {t-t )dxV' e- x ^ t ^ iC HCiV r e^^V' p T {t 

Jo 

(8) 

where 

£oPT = tz[Ho,pt], and C\pr = — [Hi , pr] • (9) 
in in 

Now, for the interaction, we adopt the form 

Hi = XY, (10) 

where X is an operator of the system and Y is an operator of the thermal bath. It is straight- 
forward to extend eq. (|10p into a more general form consisting of the sum of interactions 
(EiXiYi). 

By the standard way, we derive an explicit form of the equation as 

dp s 1 



at l h [Hs > ps] 



+ (Ji\ J* h (xe- luHs/h Xp s (t - u)e iuHs/h ^(u) - Xe- iuHs/h p s (t - u)Xe iuHs/h ^(-u) 

-e- iuHs/h X Ps {t - u)e iuHs/h X^(u) + e- iuHs/h p s (t - u)Xe iuHs/h X^{-u^j du 

+ATrBiAe ( ' _ ' o)i£o PVT(to)+A 2 TrB^ie ( ^' o)4£o f (t-to)dxV'e- x(t - to)iCo iCiV'e xit - to)lC °V' p T {t ). 

Jo 

(11) 

Here, 

X(t) = e iM s /h Xe -itH s /h_ ( 12 ) 

Because the second term on the right-hand side has the coefficient A 2 , we used the replace- 
ment 27 

ps(t) = e -(*-^sM pg ( r ) e (i-r)iHsA ( 13 ) 

Then, we have the following time-convolutionless form: 

^ = L^( P s)^-AHs,ps] 
dt in 
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+ Tr B jT ^ (XX(-u)ps(t)*(u) - X Ps (t)X(-u)^(-u)) 

-X(-u) PS (t)X^(u) + ps(t)X(-u)XV(-u)) du 

+ATr B iAe ( ^' o)i£o PVT(to)+A 2 TVBi/:ie ( ^' o)i£o [ (t-1v)dj&e-< t -*^iCiP , ^*^7*pTfai) 

J o 

(14) 

^(i) is the autocorrelation function of V in the thermal bath expressed as 

V(t) = TT B e Lm * lri Ye- im * lh Ye- m * /Z B = (Y(t)Y) B , (15) 

where 

Z B = Tr B e-^ B , (16) 

and (• • • )b denotes the average in the equilibrium state of the thermal bath. 

As natural properties of the thermal bath, we assume the following properties of the 
correlation 

lim #(i) = 0, (17) 
and the Kubo-Martin-Schwinger relation 

m(t) = ^(-t-ih(3). (18) 
Because the equilibrium state of the total system expressed as 

p ^ = L__ : z = Tie -pH? ( 19) 

is the steady state of the time evolution 

d -^T = -An^e-^/z\ =0, (20) 
at in L J 

the reduced density matrix of the system 

pf ee Tr B /4 q . (21) 

must be the steady state of the time evolution (|14j) . In the following, we study the contributions 
to the time evolution of p^ from the three parts up to the second order of A, and explicitly 
study how eq,([2U|) holds. 

Here, let us briefly study the relation between the present way of discussion and the time- 
convolutionless (TCL) formalism. In order to study non-Markov noise, the TCL formalism 
has been introduced. 26,27 The time-convolution term in eq. ([8]) was transformed into the TCL 
form 



dt 

where 



& pf{t) = tVC P f{t) - mpfit) + iVCe(t)e lV ' a V'p T (0), (22) 



V'{S(t)R(-t)-l} N 

^-^^TT^spHR}/ (23) 
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R(t) = exp fr 



S(t) = exp_ 



dre 



-iCo 



Lie 



dTV'e iCoT L ie 



(24) 



(25) 



Here, exp^_ and exp^ denote the time-ordered exponential functions. Up to the second order 
of A, this formula gives eq. (|14p . which has the TCL form. Thus, the argument in the present 
paper is for the TCL form. If we study the time-convolution form eq. (11 If) . we have a difference 
on the order of A 4 . Up to the order of A 2 , they have the same form. Therefore, the present 
argument up to the order of A 2 is good for both cases. However, as we will discuss in the next 
section, the solution of the quantum master equation has ambiguity on the order of A 2 . Thus, 
the choice of higher-order terms gives different solutions. This point will be studied in the 
next section. 

2.2 Reduced density matrix for the equilibrium state 

First, we obtain the reduced density matrix of the system for the equilibrium state up to 

the second order of A. We expand as 

e -/3(H s +H B ) 



eq 



1-0 dxe x ^ +n ^Hie- x ^ +H ^ 



+P 2 f 1 dx rdye^ +WB )%-^ +HB )e»^ +WB )%^+ WB ) + 0(A 3 
Jo Jo 



where Z' = Z$Z-q — Xa + A 2 6 with Z$ = Tr$e , Z-q = Tree 



and 



a = Tr/3 C dxe xms+nB) nie~ x ^ Hs+HB) , 
Jo 

Tr/3 2 C dx [\y e ^s+n B ) Hie - x ms+n B ) e ym s +n B ) nie -ym s +n B )^ 
Jo Jo 



(26) 



(27) 



(28) 



where Trs is the trace operation of the degree of freedom of the system. We may include 
the average of the interaction r YiH\e~^ B /Z-q into the system Hamiltonian, i.e., Tis — ► Tls + 
TtbHic-^/Zb, and then we can generally assume that 

TvHie~ mB = 0. (29) 

In this case, 

a = 0, (30) 



and thus, up to the second order of A, the is expressed in the following form: 

-/3{H S +H B ) / fl 

1 _ x 2 b - 0X / dxe x ^ +H ^Hie- x ^ +n ^ 



eq 

Pt 



ZsZb 



+/? 2 A 2 £ dx £ d ye a; ' 3 ^+WB) Wie -^(Hs+WB) e y/3(Hs+HB) Wie -y/3(W s+ W B )^ + 0(a3)) (31) 
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2.3 Steady state 

Now, we substitute eq. (|3ip into eq. (|14p. The terms on the order of 0(A 2 ) from the first, 
second, and third terms on the right-hand side of eq. (|14|) are given by 

g X 2 -/3Hs rl fx 

Ai = tL- — dx dy [n s ,X(-ih(3x)X(-ihpy)} * (iH0(x - y)) , (32) 

z s Jo JO 



) ° (xX(-u)e- /3Ws *(«)-Xe-^ Ks X(-n)*(-u) 

-X(-«)e _/3Ws X*(w) + e~ /3Ws X(-n)X^(-n)) (33) 



and 



—8X C r i 
A 3 = — — / dx X,e _/3Ws X(-t + t -ify3x) ¥(-t + «o - iTi/fa), (34) 
iftZs 7 L J 

respectively (see Appendix A). 

As we have proven in Appendix B, we find 

Ai + A 2 + A 3 = 0. (35) 

Namely, up to the second order of A, the equilibrium density matrix of the system p^?^ 
satisfies 

Cp c ^ 2) = 0. (36) 

Thus, we explicitly confirmed that the steady state of the equation of motion, eq. (I14j) . is 
given by the equilibrium state of the system under the influence of the interaction with the 
thermal bath, and that it differs from the equilibrium state of the system alone: 

fiP ^ e -W S/Zs . (37) 

It should be noted that the third term gives a non zero contribution A3 to maintain the 
modified steady state, eq. (|36|) . Often, the initial density matrix is given by a product of the 
form px(O) = Ps(O)PB) where the contribution of the third term is zero. In such cases, eq. 
(|35p is not satisfied. This means that the steady state of the total system is not given by 
a product state. Whenever we obtain a reduced density matrix for the steady state of the 
system, information on the total density matrix is lost, and then eq. (|14p holds no more. This 
apparent discrepancy is attributed to the fact that eq. (|14[) does not have a closed form for 
ps- As far as we are concerned with the reduced density matrix, we may say that the density 
matrix relaxes to pg q , but this does not mean the total density matrix relaxes to px = PsPb- 

We assumed that the autocorrelation function in the bath decays to zero in a long time 
limit, i.e., eq. (|17p . Thus, in a long time compared to the relaxation time of the autocorrelation 
function, we may ignore the third term. Then, eq. ()14[) apparently has a closed form, and we 
may regard it as an equation of motion of the reduced density matrix. Usually, we derive a 
quantum master equation in this limit. However, we must be careful with the condition to 
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ignore the third term. For example, the master equation composed only of the first and second 
terms cannot be applied to the time evolution within the relaxation time of the thermal bath. 

3. Master Equation for the Equilibrium 

Now, we consider relation between conventional quantum master equations and eq. (|14p . 
When we study a relaxation process for the equilibrium of the system (S), we usually use a 
kind of quantum master equation, which is a closed equation of motion of the density matrix 
of the system. In eq. (I14p . the third term disappears in the limit to ~°° because of eq. 
(|17p . We can regard the remaining equation as an equation of motion of p$. In this way, we 
can derive a quantum master equation. For example, we used the following compact form of 
the master equation to study several time dependences of the quantum state. 35 

^ = 4 2) ( PS ) = ^ [n s , Ps ] - (xRps(t) - R Ps (t)X - Xps(t)Ri + p s (t)rfx) . (38) 

Here, the matrix R is defined by the matrix elements given by 



/ „ „ -> T I El— E m \ J ( Em— El 

Rim ~ - % X lm * I — jj— ] " l X lm e p (Et -Em) _ 1 ' ( 39 ) 







where we use the basis that diagonalizes the system Hamiltonian Tis 

H s \m) = E m \m). (40) 

I{oS) is given by the form 

I(uj) = 9(u)I (u) (41) 

with a function Iq(u). For the arbitrary autocorrelation function Vl/(i), which satisfies the 
Kubo-Martin-Schwinger relation (|18p . we can assign eq. (|39p . 

It is easy to verify that the equilibrium of the system gives the steady state of this equation: 

4Vs q °) = 0. (42) 

In the steady state of the above master equation, no effects of the interaction with the thermal 
bath do appear, i.e., eq. (f4*2|) . although they exist in eq. (j36|) . This property (|4*2|) originates 
from the fact that we neglected the contribution of the principal value integral coming from 
the relation 



oo ■ 



e lut dt = ir5(v)+P-. (43) 



in the evaluation of the second term of eq. (|14p when we derived eq. (|38p . 35 

Now let us study the contribution of the principal value integral. In the previous section, 
we have shown that p^ 2 ^ , which includes the effects of the interaction with the thermal bath, 
is the steady state of the time evolution eq. (|14p . For it, the three contributions canceled out 
Ai + A2 + A3 = and eq. ([55]) held. Taking into account the property (|17|). the third term 
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disappears (A3 = 0) in the limit to — > —00. Thus, in this limit, we study a quantum master 
equation given by 

A x 2 



+ ( -i y o (xx(-u) Ps m(u) - x PS (t)x(-u)v(-u)) 

-X(-u) Ps (t)X^(u) + p s (t)X(-u)X^(-u)) du. (44) 



Here, we have 



Ai + A 2 = 0. (45) 

The contribution of the second term A2 consists of two parts: one from the term tt8(v), which 
we call "the real part" (A 2 r), and one from the principal value integral, which we call "the 
imaginary part" (A21): 

A 2 = A 2R + A 2 i. (46) 

In Appendix C, we show that the contribution of the principal value integral A 2 i is the same 
as — Ai: 

Ai = -A21 (47) 

up to the order A 2 . Because the equilibrium density matrix of the system itself commutes with 
H.S, the extra dynamics of the first term of the time evolution equation flHJ), i.e., Ai, is due 
to the interaction with the thermal bath. This motion is canceled by A 2 j. 

Thus, we conclude that the imaginary part gives the renormalization of Hamiltonian of 
the system. This effect has been pointed out as the Lamb shift in the driven damped oscillator 
system. 18 We can take into account the renormalization of the system Hamiltonian by adopting 
the contribution of the principal value integral. 

Now, we consider types of master equations. We often need a quantum master equation 
that causes the relaxation to the equilibrium of the system. In such cases, we can have a 
desirable equation by discarding the principal value part. On the other hand, if we need to 
discuss the effect of the contact with the bath, we have to choose the master equation including 
the imaginary part. We have to choose the types according to the purpose of the study. 

It should be noted that even if the equation of motion is correct up to the order of A 2 , the 
diagonal part of the solution has ambiguity on the order of A 2 . Let p s t be the solution of the 
master equation (|H|) 

4g(^) = 0. (48) 
If we add the arbitrary traceless diagonal operator W to p s t in the form 

pi = p s t + X 2 W. TrW = 0, (49) 
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eq. (I44p still holds up to the order of A 2 , because for any traceless diagonal operator W 

A 1 = ^[H S ,W] = 0, (50) 

and the contribution of X 2 W to A2 is on the order of A 4 . Therefore, although p^ 1 ^ satisfies 
eq. dlH) up to the order of A 2 , it does not necessarily agree with p st on the order of A 2 . 
The diagonal elements of p s t and p^ 1 ^ coincide in the 0th order of A, and their off-diagonal 
elements agree with each other on the order of A 2 , which is the leading order. The part of 
order of A 2 of the diagonal elements is related to the 0(A 2 ) terms of off-diagonal elements. 
Therefore, the 0(A 2 ) term of the steady solution of the equation depends on the 0(A ) term 
of the equation. The choice of TC (eq. (fTTj) ) or TCL (eq. (fT4l) ) gives different 0(A 4 ) terms, 
and thus they give different solutions on the order of A 2 . We may compare L^q(p^ 2 ^) with 
that obtained in the TC formalism (i.e., from eq. (jlip ). However, the evaluation of the TC 
formalism requires integration and thus it is difficult to estimate it explicitly. 

Finally, let us consider the relationship between eqs. (|14h and (138 j) . Although the reduced 
density matrices such as p^° and Pg q ^ are the solutions of the master equation, any decoupled 
initial condition in the form p^° pb cannot satisfy eq. (|14p as we have discussed in the previous 
section. Thus, it is a difficult problem to relate the master equation (|38|) and the original 
relation, eq. (|14p . Here, we may say that in the master equation we adopt the same mechanism 
of relaxation as that in eq. (fl4"|) . but they are essentially different types of equations. 

4. Summary and Discussion 

We have studied the dynamics of the coarse-grained density matrix p$. The interaction 
between the system and the thermal bath causes a modification from the equilibrium density 
matrix of the system itself. We explicitly calculated the roles of the three terms of the equa- 
tion of motion ()14[) in the steady state. We pointed out that the third term expressing the 
contribution of the initial condition plays an important role, and that we cannot ignore it in 
the study of the stationary state. 

It is often mentioned that the diagonal element of £(pg q0 ) is zero up to the second order 
of A. Although this is true, it does not help in the study of the steady (equilibrium) state of 
the time evolution. As time goes, the time evolution given by eq. (|14p always brings the state 
to p^ 2 \ but not to /9 s q0 (= e _/3s /Zs). We also pointed out that eq. (fT4|) . or more generally 
eq. ([6]), cannot be regarded as an equation of motion of ps, but it is part of the equations of 
motion for the total density matrix. 

The relationship between eq. (j!4[) and a conventional master equation was also discussed. 
Generally, a master equation is regarded as an equation obtained in the limit of to — > —00 
in eq. (|14p . There, the roles of the real part (delta function) and imaginary part (principal 
value integral) of eq. (|43p were clarified. The former drives the density matrix to that of the 
equilibrium of the system itself, and the latter modifies the system Hamiltonian due to the 
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contact with the bath. It should be noted that the latter effect is not due to the thermal 
average of the interaction part. This average has been taken into account for the system 
Hamiltonian by eq. (|29p . Thus, the modification due to the principal value integral is a kind 
of renormalization due to a dynamical process such as the Lamb shift. 

We also pointed out that although p^ 1 ^ satisfies eq. up to the order of A 2 , it does 
not necessarily agree with the solution of the equation p st on the order of A 2 . What we can 
say is the following: The diagonal elements of p s t and p s q ^ coincide on the Oth order of A, 
and their off-diagonal elements agree with each other on the order of A 2 , which is the leading 
order. 

Although the facts that we have derived are rather trivial, we hope that the explicit 
classification of the roles will clarify the structure of equations for the dissipative quantum 
dynamics. 
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Appendix A: 

The term of order 0(A 2 ) from the first term on the right-hand side of eq. (I14p is obtained 
by substituting the 0(A 2 ) term of eq. (f3Tj) : 

-\7is,(3 2 - — / dx / dye a;/3 ^ +w B) Wie -^s+WB) e i//3(Hs + WB) Wie -y/3(Hs+WB)i 

*n Z s Z B J J 

(A-1) 

If we substitute TL\ = XY (eq. (fTO]) ). we directly obtain Ai in eq. ([32]) . 

Next, the term of order 0(A 2 ) from the second term of eq. (|3ip . which gives A2 in eq. 
([33]) . is obtained by replacing pg in the second term of eq. (JT3J) with e~^ Ws /Zg. 

Finally, we consider the contribution of the third term. The third term on the right-hand 
side of eq. © is expressed by 

ATrBi£ie('-' ( ') 7 " i£o7, ' +A (*-* ^' i/:i7, >'pT(to). (A-2) 

Using the formula 

e A+XB =e A^ 1 + x J* dx e -Ax Be Ax\ + q( A 2) ( A . 3 ) 

for any set of operators A and B, the third term (eq. (IA-2jl ) up to the second order of A is 
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given by 

ATrB^ie (t - t(l)7 " i£oP ' (l + X j\t - t )dx e -x(t-t )V>iCoV' v t i£ir r e x(t-to)VU£ V'\ p' pT (t ) +O (A 3 ). 

(A-4) 

Using the relations 

Tt b V = Tt b 

V'iCoV = V'iCo 

e (t-t )V'iC _ p' e (t-t )iC ^ ' ' 

ViCiV =ViCi (from eq. (ggp ), 

eq. (|A-4|) becomes 

ATr B «/:ie ( ^* o)i/:o ^1 + A^ (i - t )dx p' e -<t-to)iCopf iClV > e x(t-t )iCo\ p'p T (t ) + 0(A 3 ). 

(A6) 

Because of the relation 

V'iCiP' = -ViCi + %C{P\ (A-7) 

for any n, 

= (A-8) 

holds. Thus, eq. (|A-6|) becomes 

ATr B i£ie ( *-* o)i£o 7'VT(to)+A 2 TrBiAe ( '"' o)i£o / 1 (t-t )^'e" :c(t ~ to)i£o i/:i^ / e x( *~* o)i£o PVT(to)- 

(A9) 

Now, we substitute /Ot(^o) of the form of eq. (f3Tj) into eq. (|A-9p . In the second term, we take 
only the 0-th order of A, and then the second term disappears because 

V ' e -/3{Hs+H B ) = Q ( A . 1Q ) 

In the first term, the term of the first order of A contributes to the 0(A 2 ) term as 

- ^!^Tr B iAe (t - t())i£o 7?V^ s+WB ) C dx ^3+^)^-^+^). (A . n) 
ZsZb Jo 

Since Ve'^Hi = 0, 

- ^LL^ic^-t^e-^^ f 1 dx e -ms+n B ) Hie ~xi3(ns+n B )_ (A . 12) 
ZsZb Jo 

We substitute Hi = XY, and then eq. (|A-12jl becomes 
\ 2 (3 



ihZ s Z B 



Tr B [AF, e - /3(Ws+HB) f dxX(-t + to-ih(ix)Y{-t + to-ih(3x) 
Jo 



\2a f -l -{3Hs fe-PH* 

1 dx [X, ——X(-t + t - ihPx)]Tr B — — Y(-t + t - ihf3x)Y ) . (A-13) 



ih Jo 1 z s V z B 

Noting that 



e -m B 



Tr B — — Y(-t + t - ih(3x)Y = ^(-t + t - ih/3x), (A-14) 
Zb 



eq. (IA-13j) corresponds to eq. eq. (([34 
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Appendix B: 

From now on, we derive 

L^(p c s q(2) ) =A! + A 2 + A 3 = 0. (B-l) 
In order to prove eq. (|B-l|h first we prove 

(Ai + A 2 + A 3 ) t=t0 = 0, (B-2) 

and 

^(A 1 + A 2 + A 3 ) = 0. (B-3) 

At t = t , 

A 2 = 0, (B-4) 

A 3 = ^— T T B [n h e-^ ns+n ^ f 1 dxe x ^+ n ^H I e- x ^+' H ^]. (B-5) 

ihZ s Z B Jo 

Here, it should be noted that Ai is time- independent and A 3 is equal to eq. (|A-12|) . Now we 
reform A 3 using the relation 

[A, e B f dx e~ Bx Ae Bx } + [B, e B f dx f dy e~ Bx Ae Bx e~ By Ae By ] = 0, (B-6) 
Jo Jo Jo 

which is derived from the second order of A in the expansion of 

[XA + B,e XA+B }=0, (B-7) 

and we have 

A 3 (t = t ) = -^^Tr B [Ws, 

Jo Jo 

= -Ai. (B-8) 

Therefore, eq. (jF3-2[) is proved. 

Next, with a straightforward calculation, we find that 

d 3\ 2 f 1 d ( e"^ 5 1 

— A 3 = -Z— I dx — {[X,——X(-t + t Q -ihPx)]y(-t + t -ihl3x)\ 
dt in Jo dt { Zs J 

BX 2 f 1 1 d f e _/3Ws 1 
= -^7z- / dx——\[X, ——X{-t + t - ihPx)]*(-t + t - iH0x) } 
in Jo inp ax l Zs ) 

= - y s { x <- ( + «e""T(< - *o) - [X, e-' m "X(-t + to)M-t + to)} 

- -I A2 ' (B ' 9) 

where we used the KMS relation (fTHj) . Thus, we proved eq. (|B-3D . 

Until now, we assumed the form Hi = XY, but it is not necessary to restrict this form of 
Hi to show eq. (iBdjl . 
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Appendix C: 

In this appendix, we prove that the principal value integral gives the A 2 term in the steady 
state of the time evolution equation (|14p when to = — oo, which is the case of the quantum 
master equation in the reference. 35 We put 



$(u) = — I dt^(t) e - iujt 
2vr J_oo 



(CI) 

in the second term of eq. (114j) and use eq. (143 p . The contribution of the principal value integral 
is 



\2 rco 

Lm 



rnPmn (t) H X k ipi m (t)X r , 



U) — LU r , 



where 



n (AY X *M X n <t)X 



CO - UJfr 



Ei — E m 



If we substitute p s (t) = e'^/Zs, eq. (RT2|) becomes 

/oo 
-OO 7 



^2 /■oo 



U - UOlr 



(C2) 



(C3) 



x X r-W 



UJ - UJlk 

From eq. (|C-ip . we obtain 



1 



3C 



Since 



-oo w-o;; n 27T J _ 00 



oo g—iujt 



oo g—iujt 



duj 



ctttf(t)P / dw 



-ivr (i > 0) 
Z7T (i < 0) 



-oo w _ w in 



If 1 f 00 

— / df V(t)e- Wlnt ■ m + — / dt^(t)e^™'-(- 
2vr y_oo 27T J 



(C-4) 



(C5) 



(C6) 



-27T 



(C7) 



Therefore, the first term of eq. (|C-4p is 

i2 /"oo 



A 2 



2ft 2 Z s 

If we assume eq. (fTT|) . we obtain 



/■oo 

/ dt V(*(t)e -<w, »* -^(t-ih[3)e iuJlnt )X kl X ln e-P En . (C- 
Jo { 



dtV(t-ih(3)e 



iuilnt 



dt^{t)e iullnt e- phwin +ih(3 



1 I dxV(- 
Jo 



-ih(3x)e 



pfajinX-PTiuJir, 



(C9) 
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by changing the path of integral. Therefore, 

\2 roc . . 

A x = -— ^ / dtY J {^{t)e- i ^ t X kl X ln e-^ - ^e^'X^e^ j 
s Jo l 

+ % Mr I dx J2 *{-ih(3x)e^- E ^X kl X ln e-P E K (C-10) 



i 

Similarly, the second term of eq. (|CJ-4[) is expressed as 

\2 roc , . 

M = ~WZ~ S J * S {*(t)e- iuJnlt X kl X ln e-^ - ^e^^X^X^e-^] 

+ ^ fdt2*(-i^)^"^H4e-^. (Oil) 
Therefore, the contribution of the first and second terms of eq. (|C-4p is 

Ai + A 2 = *g£ dx Y,i ^>{-ih(ix)X kl e-^e^ x X ln e-^ 

= iff to dx {Xe-^X(-ihPx)) kn *(-ih/3x). (C-12) 
Similarly, the contributions of the third and fourth terms of eq. (|C-4|) are 

A 3 + A A = - 1 -^ [ dx(e- (SHs X(-ihPx)x) V(-ihpx). (C-13) 
nZs Jo v 

After all, the contribution of the principal value integral is 

7/3 a 2 r 1 

^ = ^1 + ^2 + ^3 + ^4 = -^/ dx[X,e _/3Ws X(-i^x)]^(-i^x), (C-14) 

Jo 

which is the same as As(i = to)- 

Because of eqs. ()B-3p and (|B-4[) . we found that A 3 (i = t ) = —^x- Taking into account 
that Aj does not depend on ig, we conclude that the contribution of the principal value 
integral A in the limit to — * — oo is the same as — Ai. 
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